The scale-free texture of the fast solar wind. 
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The higher order statistics of magnetic field magnitude fluctuations in the fast quiet solar wind 
are quantified systematically scale-by-scale for the first time. We find a single global non-Gaussian 
scale free behaviour from minutes to over 5 hours. This spans the signature of an inertial range of 
magnetohydrodynamic (MHD) turbulence and a ~1// range in magnetic field components. This 
global scaling in field magnitude fluctuations is an intrinsic component of the underlying texture of 
the solar wind and it suggests a single stochastic process for magnetic field magnitude fluctuations 
operating across the full range of MHD time scales supported by the solar wind. Intriguingly, the 
magnetic field and velocity components show scale dependent dynamic alignment outside of the 
inertial range. 
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In-situ spacecraft observations of plasma parameters 
in the interplanetary high Reynolds number solar wind 
flow [l| provide time series over several decades that are 
ideally suited to studies of turbulence and of other phe- 
nomena that generate statistical scaling. These observa- 
tions show a power law range of power spectral density 
(PSD) in magnetic field and velocity components which 
is associated with an inertial range of magnetohydrody- 
namic (MHD) turbulence (see e.g.Q)- This inertial range 
phenomenology is seen from ion kinetic scales up to an 
outer scale at which there is a cross-over to a ~1 / / power 
law range of the PSD in the field components of coronal 
origin [l|, [|[ . 

There is debate as to which aspects of these obser- 
vations should be attributed to MHD turbulence. The 
dominant velocity and magnetic field components fluctu- 
ations are found in the inertial range of scales and these 
do exhibit statistical signatures consistent with evolving 
MHD turbulence: intermittency [H, non-Gaussian statis- 
tics § and anisotropy [f| 0] ■ This inertial range extends 
to lower frequencies both with increasing solar distance 
and in slow as compared to fast solar wind suggesting 
an actively evolving turbulent cascade @, Q. However, 
such a cascade evolves in the presence of fluctuations of 
coronal origin, the 'texture' of the solar wind [l(| fllj . 
This non-trivial texture, in the absence of turbulence, 
can be inferred, at least on longer scales, from the ob- 
served ~ 1// fluctuations in the magnetic field compo- 
nents. These show latitudinal variation and statistical 
scaling distinct from that of the inertial range 0, [H, . 

Since the earliest in-situ observations [l4| the relative 
importance of turbulence in these signatures has been 
debated. Recently, it has been argued that the range of 
PSD scaling exponents routinely measured in the iner- 
tial range for magnetic field components can arise from a 
non-evolving set of discontinuities [l5| . However, distinct 
physical mechanisms can share the same power law ex- 
ponent in the PSD |16| ■ In the context of Alfvenic MHD 
turbulence in the solar wind, where discussion naturally 



centres on the components of velocity and magnetic field, 
such discontinuities can arise from turbulence or be part 
of the solar wind texture. There is no unique procedure to 
distinguish these two physical signatures. To understand 
the texture of the solar wind it is thus informative to 
consider scalar parameters alongside vector components. 

Intriguingly, in the solar wind, scalar plasma parame- 
ters such as magnetic field magnitude |B|, energy density 
B and density n show scaling and this is also a feature 
of the texture of the solar wind. A power law range in 
the PSD of |B| is seen to extend through the inertial 
range into the ~ 1// range of frequencies but at 1 AU 
is typically a decade lower in power compared to that of 
the magnetic field components [l7| . Scalar magnetic 
field and plasma parameters can also show scaling in long 
time bulk hourly averages [lH and aggregate statistics up 
to tens of hours [ljj, l20(. Pioneering work with HELIOS 
[2ll I22I [23[ showed that |B| power spectral exponent is 
roughly invariant with radial distance at frequencies be- 
low 10~ 2 Hz; a 'flattening' at higher frequencies seen in 
the inner heliosphere is not seen at 1AU (however our 
analysis will not approach these high frequencies). There 
is an admixture of compressive fluctuations and pressure 
balanced structures [24[; |B| and n do not simply advect 
together as passive scalars |25| . 

Distinct physical mechanisms that generate sc aling can 
share the same power law exponent in the PSD |16| and 
thus statistical scaling of the higher order moments is an 
essential tool needed to distinguish them. In this Let- 
ter we perform the first systematic scale-by-scale study 
of the statistical scaling properties of |B| fluctuations 
in extended intervals of fast, quiet solar wind. We find 
global scale-free behaviour in |B fluctuations which oc- 
curs alongside, but is quite distinct from, that of the 
components of magnetic field, and extends through both 
the inertial and ~1// ranges of temporal scales. Thus we 
establish scaling is an intrinsic feature of the texture of 
the solar wind operating through the inertial and ~1// 
ranges. This suggests that a single stochastic process 
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for magnetic field magnitude fluctuations is operating or 
has operated across this full range of MHD time scales 
supported by the solar wind. 

A corollary of this result is that it provides a natural 
laboratory to test proposed measures of in-situ MHD tur- 
bulence, since within the same dataset fluctuations may 
be generated by turbulence (the inertial range) and on 
longer scales may be of solar origin. Scaling exponents 
predicted by theories of turbulence are difficult to deter- 
mine accurately in data [27l . [28| , hence the attraction to 
test signatures such as scale dependent dynamic align- 
ment [29| . Direct numerical simulations (DNS) [3(1 l3l| 
have suggested such alignment [32j . in which turbulent 
fluctuations of velocity and magnetic field progressively 
align in the cascade, to be a signature of anisotropic 
MHD turbulence. This has recently attracted contro- 
versy [HI, HH, which we will address directly from the 
observations. 

We use in-situ observations of fast solar wind from 
the ACE and Ulysses spacecraft (data obtained from the 
CDAWeb site http://cdaweb.gsfc.nasa.gov/). ACE is at 
~ 1 AU in the ecliptic, whereas for the intervals under 
study, Ulysses was at a heliospheric latitude above 70° 
at ~ 2.2 AU; thus we can compare solar wind of distinct 
coronal origin. We selected continuous intervals of quiet 
fast solar wind flow at solar minimum which did not con- 
tain large coherent structures or large scale, long time 
secular change in the ion plasma pressure (estimated as 
nT). Fast wind intervals have been identified using the 
following criteria: (v sw ) > 550 km/s and its standard de- 
viation a(v sw ) < 50 km/s, where v sw is the speed of the 
solar wind. Five stationary fast solar wind intervals were 
identified in ACE 2007 and 2008 datasets of 64 second 
average calibrated magnetic field magnitude and plasma 
observations. The Ulysses observations correspond to 5 
intervals of 3 days duration (July 1-17) of the 1995 (solar 
minimum) north polar pass. We used 1 minute aver- 
age calibrated magnetic field data for the magnetic field 
analysis only. Ulysses plasma 4 minutes averaged obser- 
vations were combined with magnetic field data by se- 
lecting the 1 second cadence magnetic field records that 
match plasma observations times to the nearest second. 

All statistical measures are first computed individually 
for each interval, these are cubic spline interpolated onto 
a common temporal or frequency grid and then an av- 
erage over all 5 intervals for each of ACE and Ulysses is 
obtained. The longest time scale is fixed by that of the 
shortest data interval so that we can access time scales 
over the range ~2 minutes to ~20 hours. Practically, 
given the effect of averaging onto 64 seconds and non- 
uniformities in the time base, we do not draw strong 
conclusions from analysis on time scales shorter than ~3 
minutes. We will indicate with error bars on our plots 
the r.m.s. variation about the sample average across the 
intervals. The averaged, smoothed PSD (Welch modified 
periodogram) of the magnetic field magnitude for both 
ACE and Ulysses intervals is shown in Figure Q] along- 
side that of one of the components: B x (ACE) and Br 
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FIG. 1: Averaged PSD of magnetic field GSE component B x 
from ACE (black circle), RTN component Br from Ulysses 
(blue circle) and magnetic field magnitude: ACE-black dia- 
mond, Ulysses-blue diamond. PSDs were displaced vertically 
for clarity; dashed lines indicate time scales of 20 and 30 min- 
utes. 



(Ulysses) . These curves have been displaced vertically in 
the figure for clarity. Typically, the magnetic field com- 
ponents are a factor of 5 higher in power compared to 
|B| for these intervals and on this plot show a steeper 
power-law trend at higher frequencies with exponent, s, 
close to the often observed ps | 'inertial range' value. 
Figure [T] shows that at a time scale 20-30 minutes the 
components depart from w | inertial range scaling and 
cross-over to a spectrum. This time scale is indi- 

cated on all subsequent plots. The PSD of |B| follows 
a single power-law over the entire range with exponent 
w— 1.4, a behaviour that is distinct from that of the com- 
ponents. The Ulysses and ACE observations show essen- 
tially the same behaviour in the PSD: a transition from 
inertial range to ~ 1// scaling in the components (also 
seen in the other components, not shown) and alongside 
this, a single range of scaling in the magnitude with ex- 
ponent distinct from that of the inertial range. We now 
consider the statistical scaling in the higher moments of 
the fluctuations. On a given temporal scale r, the fluc- 
tuations are 8x{t,r) = x(t + r) — x(t) in the time se- 
ries of a given (scalar or vector) quantity x(t). In Fig- 
ure [5] we plot, for the same quantities shown in Figure 
[TJ the Kurtosis of the fluctuations K — fi^/ where 
/ife =< (Sx(t,T)— < 6x(t,r) >) k > are their central 
moments. For a Gaussian PDF K = 3 and this is in- 
dicated with a dashed horizontal line. The components 
(within variance of the average) show a smooth and scale 
free drop in Kurtosis within the inertial range, on scales 
shorter than ~ 20 — 30 minutes. At ~ 20 — 30 minutes 
there is a 'kink' at the end if the inertial range, beyond 
which the PDFs are weakly non-Gaussian (K > 3) and 
progressively approach Gaussian at scales ~ 120 minutes. 
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FIG. 2: Kurtosis of magnetic field component B x from ACE 
(black circle), component Br from Ulysses (blue circle) and 
magnetic field magnitude: ACE-black diamond, Ulysses-blue 
diamond. 
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FIG. 3: Average PSD of magnetic field magnitude (ACE- 
green triangles, Ulysses-black diamond) and density (ACE - 
red circle, Ulysses-blue square). 



The behaviour of <5|B| is clearly distinct from that of the 
components, and is far from Gaussian until much longer 
time scales, beyond ~5 hours. On scales shorter than 
~5 hours the behaviour in the Kurtosis seen in fast quiet 
solar wind streams in the ecliptic (ACE) roughly tracks 
that seen in polar outflows (Ulysses); fluctuations of |B| 
are not sensitive to the distinct coronal origin of these 
flows until the longest temporal scales are reached. We 
note that, in principle, more information can be obtained 
from signed structure functions of odd moments, how- 
ever, consistent with [H, [36J, we find that these signed 
quantities do not produce homogeneous results across all 
intervals studied here, therefore it is not appropriate to 
average these. 

Figures Q] and [2] together establish a quite remarkable 
result: there is a single scale-free behaviour in fluctua- 
tions of |B| spanning three decades in temporal scales 
from ~ 2 minutes to over 5 hours. The scale-free signa- 
tures extend through both the inertial and ranges 
of temporal scales seen in these same intervals in the field 
components. This suggests a physical mechanism for |B| 
scaling across all these scales that is distinct from the in- 
situ Aflvenic MHD turbulence that is driving scaling in 
the components. Compressive fluctuations and pressure 
balanced structures [24| are natural candidates for mag- 
netic field magnitude fluctuations-they can in principle 
be dynamically evolving or passively advecting, having 
been generated in the formation of the solar wind. In 
Figure [3] we directly compare the PSD of |B| from ACE 
and Ulysses as shown in Figure [1] with that of the den- 
sity. The scaling of the density closely follows that of |B| 
in the ~1// range, but departs on time scales shorter 
than 20-30 minutes, that is, where the inertial range is 
seen in the magnetic field components. In the inertial 
range, the power in the density fluctuations is enhanced 



above that which would arise from the scaling seen on 
longer time scales. Essentially the same behaviour is 
seen in both ACE and Ulysses and thus again is not 
strongly sensitive to the region of coronal origin. The 
single range of scaling in |B| fluctuations that extends 
through both the inertial and ranges of scale thus 

does not straightforwardly correspond to the behaviour 
of the density fluctuations. Thus one cannot invoke a sin- 
gle process to generate the scaling of fluctuations in both 
|B| and density unless some physical process operating 
in the inertial range acts to enhance density fluctuations. 

Finally, we will look more closely at the behaviour 
of the components of fluctuations across both the iner- 
tial and ~ 1// ranges of scale. We explore this idea by 
performing a dynamic alignment analysis across the full 
range of scales available in our chosen intervals. Dynamic 
alignment, calculated scale-by-scale as the relative angle 
between the vector fluctuations Sv± and Sb± [2{| [32j , 
has been proposed as a signature of the turbulent cas- 
cade [32j |. Perpendicular directions are taken with re- 
spect to the mean local magnetic field, which at time 
t is given by an average (B(f)) computed over a scale 
dependent interval [t — r, t + r] . Perpendicular fluctu- 
ations are then obtained from the following expressions: 
5v± = <5v — (dv • b)b and <5b^ = <5b (<5b • b)b, where the 
unit vector b = (5B/B. The average angle of alignment 
is: 



0(t) = arcsin 



dv± A Sbj 

(Sv±) (6b ± 



(1) 



In Figure [4] we plot this average angle O(t) as a func- 
tion of scale t for the same intervals shown in Figures 
[T][3J On such a plot, progressively aligning field and flow 
fluctuations should show a scale free decrease in (O) with 
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FIG. 4: Average angle of perpendicular fluctuations in ve- 
locity and magnetic field ACE-black diamond, Ulysses-blue 
square. An asymptote of ~ r - 25 jg indicated on the plot. 



decreasing scale r within the inertial range; as the cas- 
cade progresses from large to small scales the field and 
flow progressively align 32]. However, Figure[4]show this 
behaviour most strongly on temporal scales r > 100 min- 
utes, well into the range of scales. The alignment 
flattens on shorter scales, reaching its minimum at the 
upper end of the inertial range at 20-30 minutes. Indeed, 
the variation of (0) with decreasing scale r asymptotes 
to within errors of that predicted for in-situ turbulence 



[HJ ■ Thus "dynamic alignment" , as quantified by ([TJ , is 
not a clear discriminator of turbulence in the solar wind, 
but intriguingly, in the absence of in-situ MHD turbu- 
lence, that is, in the ~1// range it is coincident with 
scaling in fluctuations of |B|. 

We have established that there is a scale- free behaviour 
in |B| fluctuations that occurs alongside, but is quite dis- 
tinct from, that of the components of magnetic field. It 
extends over three decades through both the inertial and 
~1// ranges of temporal scales and is seen in fast quiet 
solar wind of different coronal origin. The most parsi- 
monious description of how such a single scaling range in 
|B| could arise is that a single process operates, or has 
operated, over all these scales. Observational evidence of 
incompressible MHD turbulence in the solar wind must 
thus be understood in the context of the global scaling 
of the 'texture' [l^, [H| of the solar wind. We have also 
shown that scale dependent 'dynamic alignment' is not a 
clear discriminator of turbulence in the solar wind. This 
'pseudo-dynamic alignment', taken alongside the scaling 
seen in the magnitude of magnetic field, may however 
provide an insight into the generation of the flux tube 
texture of the solar wind. It may capture some physics 
of the generation of the solar wind, reflecting the manner 
in which magnetic helicity is injected via photospheric 
fields that are fractal [37| . 
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